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What flows in the chiral magnetic effect? 
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To address a question of whether the chiral magnetic current is a static polarization or a genuine 
flow of charged particles, we elucidate the numerical formulation to simulate the net production 
of right-handed particles and anomalous currents with CP-breaking background fields which cause 
an imbalance between particles and anti-particles. For a concrete demonstration we numerically 
impose pulsed electric and magnetic fields to confirm our answer to the question that the produced 
net particles flow in the dynamical chiral magnetic effect. The rate for the particle production and 
the chiral magnetic current generation is quantitatively consistent with the axial anomaly, while 
they appear with a finite response time. We emphasize the importance to quantify the response 
time that would suppress observable effects of the anomalous current. 


I. INTRODUCTION 

In many quantum problems in physics it is highly de¬ 
manded to establish a numerical framework to simulate 
full dynamical processes of particle production out of 
equilibrium. Unsettled and important problems include 
the leptogenesis and the baryogenesis for the explanation 
of the baryon asymmetry of the Universe (BAU). It is 
well-known that Sakharov’s three conditions are needed 
for the BAU; namely, 13-breaking process, C and CV vi¬ 
olation, and out-of-equilibrium. The particle production 
on top of the electroweak sphaleron [Tj accommodates a 
process with A(13 + L) ^ 0, which means that B + L 
decays toward chemical-equilibrium [and this is why the 
idea based on the SU(5) grand unified theory [2] does 
not work out for the BAU]. Therefore, if the leptogenesis 
with some (B — L)-breaking process beyond the Standard 
Model (see Ref. [3] for example) generates L / 0 initially, 
it would amount to B 0 finally in chemical-equilibrium. 
The particle production associated with the sphaleron 
transition has a character of topological invariance and 
the index theorem relates the topological number to the 
change in the particle number. We must recall, how¬ 
ever, that the sphaleron is a special static configuration 
at saddle-point and for more general and non-topological 
gauge configurations with C- and CP- odd components 
we have to consider microscopic simulations with Weyl 
fermions. Also, to compute the momentum spectrum 
of produced particles with arbitrary gauge backgrounds, 
there is no longer mathematical elegance and some brute- 
force tactics with numerical simulations is indispensable. 

The sphaleron transition rate is proportional to T 4 
where T is the temperature @], and so such topologi¬ 
cal fluctuations should be abundant also in the strong 
interaction when the physical system is heated up to 
T > Aqcd as argued in Ref. [5]. Such high tempera¬ 
ture is indeed realized in the quark-gluon plasma cre¬ 
ated in the relativistic heavy-ion collision (HIC) exper¬ 
iment. Although the strong interaction does not break 
CP (except for negligibly small strong 0), we can still 
anticipate local violation of V and CV before thermal- 


ization [BJ. This possibility is nowadays referred to com¬ 
monly as the local parity violation (LPV). Theoretical 
and experimental challenges are still ongoing about the 
LPV detection (and no direct evidence has been found 
for it). Along these lines the recognition of the inter¬ 
play with pulsed magnetic field B between two heavy 
ions was a major breakthrough 0; B on top of CV- odd 
backgrounds would induce an electric current j in par¬ 
allel to B, which is summarized in a compact formula 
of the chiral magnetic effect (CME) with chiral chemi¬ 
cal potential //5 [8]. In a realistic situation of the HIC 
we can anticipate V and CV violation not only from the 
QCD sphalerons but also from the glasma [9] which is a 
description of the initial condition of the HIC in terms 
of coherent fields. Then, without introducing /Z5 that 
has huge theoretical uncertainty (see Ref. [121 for a recent 
study with /X5), we can and should directly compute the 
anomalous currents associated with the particle produc¬ 
tion with CP-odd backgrounds composed from chromo¬ 
electric and chromo-magnetic fields un¬ 
it would be an intriguing attempt to test the CME in 
a table-top experiment in a way similar to the quantum 
Hall effect. The biggest difference between the CME and 
the Hall effect arises from the fact that the carriers of 
electric charge are not the ordinary particles but chiral 
particles; that is, a non-zero excess of the right-handed 
particle number to the left-handed particle number is re¬ 
quired m ■ It is quite non-trivial how to implement such 
a chiral battery (usually represented by /i5 7? 0) in ma¬ 
terials. Just recently the CME may have been confirmed 
in a system rather similar to the glasma, in which chiral¬ 
ity imbalance is imposed by background fields that break 
CV symmetry [Ill- 

Here, to elucidate the motivation of our present study 
clearly, let us discuss theoretical subtleties about the in¬ 
terpretation of the chiral magnetic (and separation) ef¬ 
fect. 

1. What carries the electric charge of the current? 
In the derivation using the Chern-Simons-Maxwell 
theory [14] the chiral magnetic current appears 
in Ampere’s law, but it takes a similar form as 
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Maxwell’s displacement current. We know that 
Maxwell’s displacement current is a source for the 
magnetic field but no charge carrier flows (though 
the Poynting vector shows a flow of electromagnetic 
energy) as pointed out in Ref. fI5j . 

2. What is the momentum spectrum of charged parti¬ 
cles that flow as the chiral magnetic current? In the 
derivation using the thermodynamic potential [Sj a 
finite result remains after a cancellation between in¬ 
finitely large momentum contributions p z ~ ±oo, 
but it is unlikely in any experiment that particles 
with such large momenta emerge from the vacuum. 
In the kinetic derivation |16H18] a thermal distribu¬ 
tion of particles is finally assumed to retrieve the 
exact CME coefficient, which provides us with a 
theoretically correct description of the dynamical 
current. It is, however, still needed to understand 
how such a distribution of particles arises. In gen¬ 
eral with arbitrary (non-thermal) initial conditions 
the particle distribution functions are not neces¬ 
sarily thermal and then the CME current may not 
obey the standard formula, which also motivates 
the numerical simulation. 

3. The expression (f2|j|fl) ^ 0 does not always mean a 
flowing current but it may represent a polarization 
that is a static expectation value of a vector oper¬ 
ator j. where Q) represents a certain pure state of 
the physical system. (We can generalize our discus¬ 
sions to mixed states, but the consideration with a 
pure state suffices for our purpose.) In the lattice- 
QCD simulation in Euclid spacetime as in Ref. [T9], 
fi) is the Euclidean QCD vacuum, and the current 
measured with Q) should be a polarization because 
the current is a non-equilibrium and steady dynam¬ 
ical phenomenon. Actually the chiral separation 
effect (f 2 |j A |fi) oc pB has a natural interpretation 
as a spin polarization [20] (see also Ref. [2T] for 
a recent study on this). To make a possible dif¬ 
ference explicit we point out an example found in 
the computation of the chiral magnetic conductiv¬ 
ity a x (oj,p) that represents a linear response in the 
presence of spacetime dependent B(uj,p). From 
Fig. 2 of Ref. [22] it is obvious that the static limit 
a x (co = 0 ,p —> 0 ) takes a value different from the 
dynamical limit ct x (uj —> 0 ,p = 0 ). 

4. How can the particle motion be aligned to the mag¬ 
netic field? The conventional explanation (i.e., the 
correlation between the helicity and the spin align¬ 
ment under a strong B) is based on a static picture 
corresponding to the polarization phenomenon. If 
there are free right-handed fermions, say, in the 
direction perpendicular to B, they just move on 
a circle with the Larmor radius classically. The 
spin receives a torque from the Berry’s curvature 
term, but the spin cannot be aligned to B with¬ 
out dissipation. Also we make a comment that the 
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FIG. 1. Schematic view of currents induced by B y in a CP- 
odd domain realized by parallel E z and B z . The currents 
flow in all the x, y, and .z directions; j x is the anomalous Hall 
current and j v is the CME current. 


kinetic derivation may look like a problem of one- 
particle motion, but it should be justified by the 
worldline formalism on the quantum level, and then 
the (proper) time may be given a meaning different 
from the genuine time. 

5. What is the response time for the chiral magnetic 
current to get activated? It is unlikely that the 
current suddenly starts flowing as soon as B and 
/is are turned on. (This is actually an unavoid¬ 
able problem to simulate the chiral magnetic effect 
assuming a certain p§.) The current generation 
rate has been calculated only in an idealized setup, 
and the quantitative estimate of this response time 
should be crucial for experimental detection in re¬ 
alistic and thus disturbed environments. 

To answer these question, in this work, we consider 
a special setup as illustrated in Fig. [T] which is the 
schematic view of the CME setup in the HIC and in the 
condensed matter experiment. The parallel E and B (in 
the z direction in Fig. [l]) form V- and CP-odd product 
E ■ B, and the CME current j v is induced in a direction 
perpendicular to E , which is reminiscent of the anoma¬ 
lous Hall current j x . (In this sense one can regard the 
CME as a 3D extension of the Hall effect.) In the HIC 
E z and B z are Abelian projected chromo-electric and 
chromo-magnetic fields m ■ In the sphaleron transition 
these fields are provided from the Abelian projected part 
of the sphaleron gauge configuration too. 

Then, we can address the questions in a very concrete 
way and we shall present short answers in order: 

1. The produced particles from the vacuum fluctua¬ 
tions form a current. The right-handed particles 
and anti-particles have asymmetric momentum dis¬ 
tributions due to E ■ B 7 ^ 0, where the particle 
number is conserved with the opposite excess from 
the left-handed sector (and the current is doubled). 
Therefore, the chiral magnetic effect can be real¬ 
ized in a non-equilibrium environment and is most 









3 


relevant to the early-time dynamics of the HIC in¬ 
volving the particle production. 

2. Produced particles are accelerated by E and the 
momentum should have a peak around an impulse 
given by the product of the Lorentz force times 
the duration, which is the ordinary situation in the 
Schwinger mechanism with pure E. In the same 
way we can compute the distribution function with 
B ± 0 as well as E ^ 0 that is localized in mo¬ 
mentum space and so regions with infinitely large 
momenta are completely irrelevant. 

3. The particle production process is the real-time 

phenomenon out of equilibrium. In this sense, ac¬ 
cording to our understanding based on the parti¬ 
cle production, oc p,B can be interpreted 

as a dynamical flow of the chirality, which results 
simply from a different combination of the right- 
handed and the left-handed sectors of the present 
calculation. This fact supports the realization of 
non-zero quadrupole moment of the electric charge 
suggested from the chiral magnetic wave. We note 
that if were a spin polarization or a mag¬ 

netization, there would be no movement of chirality 
in space, which is not the case in our fully dynam¬ 
ical setup. 

4. The particle motion has a peculiar angle distribu¬ 
tion or asymmetric momentum distribution (which 
will be demonstrated in Fig. |2]) reflecting the config¬ 
urations of E and B. Thus, unlike the conventional 
explanation, there is no need to bend the motion 
of classical particles nor the direction of the spin 
polarization with B only. 

5. One might think that the asymmetric momentum 
distribution occurs from the beginning of the parti¬ 
cle production and so there is no delay for the topo¬ 
logical currents to flow. The plane wave, however, 
should be deformed into the Landau wave-function 
after a finite B is switched on. The response time 
is governed by the development in the amplitude of 
the Landau wave-function. In particular the onsets 
of the particle production and the current genera¬ 
tion may look different due to different rates. 


II. THEORY OF THE PARTICLE AND 
CURRENT PRODUCTION 

In this work we focus on the right-handed sector only 
and the current (net particle) density should be doubled 
(canceled) once the left-handed sector is taken into ac¬ 
count. The right-handed Weyl fermion should satisfy the 
following equation of motion: 

(ia p d tl -ea tl A p )<f> R = 0 (1) 

with a p = (l,er). We can readily construct a complete 
set of plane-wave solutions of Eq. (|T|) for the asymptotic 


states where the interaction is turned off. In general a 
constant background A M may be coupled even without 
interaction. In the present work we set a gauge that 
makes A$ = 0 for a technical reason. We can then write 
positive-energy particle solutions as (f > R = u R (p; A)e~ lp ' x 
with 


ur (p;A) 


( v 7 !P. 4 1 + Pa A 


( 2 ) 


in a certain gauge [or a convention for the overall U(l) 
phase]. Here we defined p ±A = p z F eA and the phase 
factor is: 


p w( PA ) _ Pa + iPA m 

V(Pa) 2 + (p V a) 2 ' 

We can identify the anti-particle state from </> R = 
— ia 2 ^ R , which leads to a relation: u R (p; A) = 
ur{—P',A). In the same way we can find the 
negative-energy particle and anti-particle solutions 
with (p R /R = v R / R (p; A)e +lp ' x , which results in 
v R (p;A) = -e~ l6( - p -^u R (p;-A) and v R (p; A) = 

— e - ie (p A ) UR (—p; —A). We note that u R (p; A) and 
v R (p; A) have an energy ±|p A |, while other two, 
u R {p; A) and v R (p-,A), have an energy ±\p_ A \. 

For the problem of particle and anti-particle produc¬ 
tion we evaluate an amplitude for the transition from a 
negative-energy state (with momentum p and vector po¬ 
tential A) to a positive-energy state (with momentum q 
and vector potential A/), which we can explicitly express 
as 


Pq.p — 
Pq.p = 


3 ui{q A ')e l]qA '^ +lq x „ , 0 


d d x 




g- p (x°,x) 


v t (n „,ip i lq-A'|z°+*'?' 3: 

j 3 u yq_A>^ 0 

d x ~ - /o , | - 9-p(x , x ) 

\/2\q_A' I 


( 4 ) 


for the particles and the anti-particles, respectively. We 
here introduced new functions, g p {x) and g p {x), as so¬ 
lutions of the particle and the anti-particle equations of 
motion satisfying the following negative energy boundary 
conditions: 


9 P {x) 

g p {x) 


VRkP-A)^ p - Ax °- ip '* 

V 2 \p~a\ 

u R (p A )e i \ p A x °~ ipx 
V 2 \Pa\ 


( 5 ) 


for a; 0 around the initial time in the past. It should be 
noted that the equation of motion for the anti-particle 
is not Eq. |l]) but e is replaced with —e and a p with 
d p . Finally we can express the net particle number (i.e. 
the particle number minus the anti-particle number in 
the right-handed, sector) as well as the spatial currents 
in the following manner (after dropping the zero-point 














4 


oscillation term from J°): 

d 3 q 


J° = e 


J = e 


(2tt) 

d 3 q 


3 [/(<?) - /(<?)] , 


O0 3 LI <ia 


I Q-A' 


( 6 ) 


where the distribution functions f(q) and f(q) above are 
defined with the amplitudes integrated over all incoming 
momenta, i.e., 

~ J (2,r)3 1^9 -pI 2 ’ = J (^KpI 2 - 

These functions are naturally given an interpretation as 
the distribution functions of particles and anti-particles. 
We note that J° above is a formula equivalent to that 
as utilized in Ref. [233- We also make a comment on the 
relation to the chiral kinetic theory where the current is 
defined with x and the equation of motion solves bx in 
terms of q and an additional contribution from Berry’s 
curvature. The appearance of the canonical momenta 
q A f (shifted by A') in Eq. ^ corresponds to using bx 
on the level of the equation of motion. Still, it is highly 
non-trivial how to see the analytical relation between two 
approaches, which will be an interesting future problem. 


III. SIMULATION SETUP 


We here consider pulsed fields approximated by step 
functions for a duration T. We choose the origin of the 
time so that we start solving Eq. 0 numerically from 
t = 0. Introducing a temporal profile defined by 


e(t) = 


(-T/2 < t-t 0 < T/2) 
(\t-t 0 \ > T/2) , 


( 8 ) 


we can explicitly specify the electric and magnetic fields 
that we want to consider here as 


E z (t) ~ B z {t) ~ B v {t) oc e(t) . (9) 

This choice of the temporal profile is motivated based on 
the glasrna dynamics in which both the external magnetic 
field and the chromo-fields decay within the same time 
scale of ~ 0.1 fm/c [2HI2S]- Theoretically speaking, to 
define the produced particle number uniquely, we should 
setup the asymptotic states where interaction is switched 
off. We also mention another possibility to prescribe the 
particle number in a transient state (mostly relying on an 
adiabatic approximation) as studied in Ref. [26], though 
we do not adopt it. 

It is, however, technically non-trivial how to realize the 
external fields as Eq. ([9]) strictly. Here, we will make a 
simpler choice of the gauge potentials: 

A* = B ± e(t)z - B n e(t)y , A y = 0, 

A z = -E 0 e(t) , 1 J 


where e(t) = e(t')dt'. We note that A z {t = — oo) = 0 
and A z (t = oo) = —E 0 T ^ 0 and this is why we should 
keep p 4 in all the expressions for produced particles. 
The above vector potentials actually lead to the mag¬ 
netic fields as almost desired: 

B x {t) = 0, B*(t)=B ± e(t), B z (t)=B {l e(t) , (11) 

whereas the electric fields have some extra components: 

E x {t) = - B ± e\t)z + Bue'(t)y , E y (t) = 0 , 


E z {t) =E 0 e(t) . 


( 12 ) 


We chose the gauge (10) so that we can have E v (t) = 0 


because we are interested in j y . This also indicates that 
the current j x has not only the anomalous component 
but a finite contribution induced by E°° ^ 0. 

We note that in our present calculation we neglect the 
back-reaction and so this j x would not affect j y . Be¬ 
cause we start the simulation with the null initial condi¬ 
tion with zero particle, this approximation for the back- 
reaction is expected to hold well especially for our setup 
with pulsed fields. In principle our formulation can be 
upgraded to the classical statistical simulation contain¬ 
ing back-reaction m or the Dirac equations can be com¬ 
bined with a numerical solution of Maxwell’s equations. 
The classical evolution suffices for our present purpose 
to figure out the net particle production and associated 
CME current generation. 

Let us now go into some more technical parts. It is con¬ 
venient to keep the reflection symmetry of axes in order 
to avoid non-vanishing currents due to lattice artifact. A 
perturbation theory in terms of A 1 would give a term like 
oc z, y that should be vanishing after the spatial integra¬ 
tions. This can be non-perturbatively realized on the lat¬ 
tice by the periodicity of the link variables. Alternatively, 
if we take the spatial sites n l (that gives x l = n l a with the 
lattice spacing a) from — IV, to -t-lV* for i = x,y, z , these 
terms oc z, y trivially disappear. The spatial volume is 
thus (2N X + 1) • (2 N y + 1) • (2 N z + 1) • a 3 in the present 
simulation. The corresponding momenta are discretized 
as p l a = 2nk l /(2Ni + 1). We also note that we impose 
the anti-periodic boundary condition to avoid the singu¬ 
larity at \p A \ = 0 in Eq. ([2]) or in Eq. Q. This means 
that k l should take half integral values from — iV) + 1/2 to 
+Ni + 1/2. We emphasize here that, if the spatial volume 
is large enough, it does not matter whether the bound¬ 
ary condition is periodic or anti-periodic. This is simply 
a prescription not to hit the singularity at \p A \ = 0. One 
might think that one could insert a small regulator in the 
denominator like the ie prescription. Indeed, one could 
adopt the periodic boundary condition and could take 
an average of contributions from six neighbors, p A = ±e, 
p A = ±e, p A = ±e around \p A \ = 0. As seen in mo¬ 
mentum space, however, such an averaging procedure is 
effectively identical to imposing the anti-periodic bound¬ 
ary condition. In any case we should carefully treat the 
singularity at \p A \ = 0 because this infrared singularity 
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Particle 

Anti-particle 



FIG. 2. Distribution of the produced particles and anti¬ 
particles as functions of q v and q z (and integrated with respect 
to q x ; see the text). The lattice size is N x = N y = N z = 12 
and the magnetic fields are B y = B± = Eq/2. 


is responsible for the axial anomaly, which is manifested 
in the representation with Berry’s curvature [16Hl8| . 

If we integrate over all momenta, we cannot avoid pick¬ 
ing up all contributions from the doublers and the net 
particle production should be then absolutely zero be¬ 
cause the axial anomaly is exactly canceled according 
to the Nielsen-Ninomiya theorem f2E). It is thus cru¬ 
cial to get rid of the doubler contributions properly. To 
this end we limit our momentum integration range to 
the half Brillouin zone only, i.e., k l from —TV)/ 2 + 1/2 
to +Ni/ 2 — 1/2 (which is used in Ref. [33] too). We 
note that the axial anomaly is correctly captured from 
the singular contributions in the wave-function around 
\p A \ = 0 , though it usually appears near the ultraviolet 
cutoff in a diagrammatic approach. (We have numeri¬ 
cally checked that our results are robust against small 
shifts at the momentum edges, which is also understand¬ 
able from Fig. [2} produced particles are centered with 
small momenta.) One should not be confused with the 
situation in Euclidean, i.e., static lattice simulation in 
which the whole Dirac determinant should be evaluated 
to identify the statistical weight for gauge configurations. 
For the present purpose, in contrast, we need the infor¬ 
mation on each momentum mode, and so the computa¬ 
tion is much more time consuming than evaluating the 
determinant, but it is easier to get rid of the doublers on 
the mode-by-mode basis. 

When we change the lattice size and the lattice spacing 
a, we should change a dimensionless E 0 a 2 accordingly to 
fix E 0 in the physical unit. In this work we choose a 
reference lattice size as N x = N y = N z = 8 (i.e., 17 3 lat¬ 
tice volume) and scale physical quantities such as Eqo 1 oc 
17 2 /(2IV+1) 2 with varying N x = N y = N z = N. To sim¬ 
plify the notation let us absorb e in the field definition 
and omit a hereafter. Then, we take Eq(N = 8 ) = 0.1 
and change By and B± from 0 to E 0 . We set the dura¬ 
tion as T = a/IO/Eo (that is, T = 10 for N = 8 ). Our 
choice of t 0 is t 0 = 0.6 T and we continue solving Eq. 0 
up to t = 1.2 T with time spacing At = 0.02T using the 


2nd-order Runge-Kutta method (with which we carefully 
checked that the accuracy is good enough for our present 
simulation). 

We did not quantize the field strength unlike Ref. 35]. 
Actually we compared results with unquantized and 
quantized fields (in the unit of 2 n/N x N y for B z etc) 
and found that this would make only negligible differ¬ 
ence. Because our calculation is classical for gauge fields, 
non-quantization of the field strength is harmless for our 
numerical results, and moreover, technically speaking, we 
do not have to impose a boundary condition for the gauge 
fields unless we solve the equation of motion containing 
the derivatives of the gauge fields. 

Even if By = B± = 0, a finite Eq would induce the 
pair production of particles and anti-particles under the 
Schwinger mechanism (see Ref. [30j for a comprehensive 
review). Because we deal with Weyl fermion, the pair 
production is always possible for any Eq ^ 0 , which leads 
to a finite distribution of particles with —EqT < p z < 0 
and p x = p v = 0 (and with 0 < p z < EqT for anti¬ 
particles). If By or B l is finite, this distribution is 
smeared in momentum space as clearly observed in Fig. [2] 
in which we present J(dq x /2n)f(q) and f(dq x /2n)f(q) 
as functions of q y and q z in the case with By = B± = 
Eq/2. In addition to smearing, it is evident in Fig. [2] 
that the particle is more enhanced over the anti-particle 
because of the CB-odd background, which signals for the 
net production of right-handed particles. (For the vector 
gauge theory the total particle number is conserved as a 
result of summing right-handed and left-handed contri¬ 
butions up.) 

To close this section we shall make some remarks about 
the lattice convention. Taking the sum over dimension¬ 
less momenta k -1 corresponds to the phase-space integral 
V f d 3 p/( 27t) 3 . In our definition j3 q _ p has a mass dimen¬ 
sion of V and so f(p) has a mass dimension of V too. 
Then, J° and J become dimensionless quantities. In 
these final expressions, however, the prescription to avoid 
the doubler contributions makes the phase-space volume 
smaller by a factor N 3 /(2N + l) 3 , which should be cor¬ 
rected in the end. In this paper we would not show the 
lattice version of expressions, but for example, q± A in 
Eq. <[6j) should be understood as sin[(g i + A l )a\. 


IV. NUMERICAL RESULTS 

In what follows we normalize j° = J°/V and j = J /V 
by the following quantity: 


B 2 

U ° = +7T 3 ' T ’ 


(13) 


which is the pair production rate Bq/ 47 t 3 in the 
Schwinger mechanism multiplied by the duration T, 
which gives us a natural order-of-magnitude estimate. 
Then, j°/ n o should remain to be a reasonable number, 
which is confirmed in Fig. [3] 
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FIG. 3. Net particle density j° normalized by no in Eq. (131 
as a function of Bn and B±. The lattice size is N x = N y = 
N z = 8 ,10,12 from the bottom (solid) to the top (dotted). 




FIG. 5. Distribution of the produced particles as functions 
of q v and q z (and integrated with respect to q x ) for different 
lattice sizes IT 3 , 21 3 , and 25 3 . The magnetic fields are B\\ = 
B± = Eq/2 and the duration is 0.5T. 


means E ■ B = 0, and so, no net particle production is 
allowed (even though the pair production is possible as 
long as E 0). The current is vanishing then due to the 
lack of electric carrier. 


V. LATTICE-SIZE DEPENDENCE AND THE 
RESPONSE TIME 


FIG. 4. Net current density j v normalized by no which is to 
be interpreted as the CME current. 


Figure [3] clearly evidences for a non-zero value of the 
net particle production for E B = E$B\\ ^ 0. We 
can make it sure that the net particle production is van¬ 
ishing when B\\ = 0. This should be trivially so but a 
non-symmetric treatment of coordinates would pick up 
unphysical artifact, which is not the case in the present 
simulation. In Fig. [3] the lattice size is N x = N y = N z = 
8 ,10,12 (i.e., 17 3 , 21 3 , 25 3 lattice volumes) from the bot¬ 
tom (solid) to the top (dotted). We can conclude that 
dependence is rather mild, while it is still there. In fact j° 
shows moderate dependence on B± , which results from 
the lattice artifact. In the steady system with constant 
electromagnetic fields the anomalous particle production 
rate has only weak dependence on B± and it slightly in¬ 
creases with increasing B±. We can see a tendency that 
j° at smaller (larger) B± is pushed down (up) as A) in¬ 
creases. 

Now we are going to discuss j/no obtained in our nu¬ 
merical calculation. Here let us focus only on the chiral 
magnetic current j v and postpone discussions about the 
anomalous Hall current j°° and Ohm’s current j z to a sep¬ 
arate publication M- It is interesting that our results 
in Fig. [4] surely represents a non-zero CME current as 
theoretically expected. We see that j y = 0 when either 
B || or Bj_ is vanishing. Because j v is the CME current, 
it is naturally vanishing if B± = 0. Furthermore, B || = 0 


The lattice-size dependence that we have seen is the 
ultraviolet extrapolation, for which the box size is fixed 
and the lattice spacing is decreased. Therefore, E 0 in 
the unit of the lattice spacing was rescaled to fix the ob¬ 
servables in the physical unit. In the present case the 
numerical results should be stable for the ultraviolet ex¬ 
trapolation. We can understand this from Fig. [5] As 
mentioned before, the distribution function is well local¬ 
ized in momentum space, and so the contributions from 
the high-momentum region simply attach a tail in the 
particle distribution as shown in Fig. [5j 

In our definition f(q) does not go to zero for \q\ —> oc 
when E 7^ 0 and B / 0. If the subtraction of the zero- 
point oscillation energy is normalized at q\ —> 00 we could 
remove an irrelevant shift in the distribution function. In 
principle, however, such a constant shift does not affect 
j° and j: the particle and the anti-particle contributions 
cancel for <7° and there is no contribution to the angle 
integration of J d 3 p for j. From this consideration we 
can expect better ultraviolet convergence for j because a 
constant shift vanishes for each of the particle the anti¬ 
particle parts, while j° needs a delicate cancellation. We 
will see that this expectation is indeed the case in our 
numerical simulation. 

It would be of pragmatic importance to consider j° and 
j v as a function of the pulse duration, that is, the life time 
of the external fields. So far, we fixed it as T = -y/10/Eo 
and let us vary it now. Theoretically speaking, for a 
sufficiently large duration, we should expect linear de- 
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FIG. 6. Profile of j° and j y as a function of the pulse duration 
normalized by no- The lattice size is N = 8 (solid line), 10 
(dashed line), 12 (dash-dotted line). The magnetic fields are 
chosen as By = Bj_ = Eq/2. The dotted line represents the 
analytically estimated rates for the particle production and 
the current generation (with a common offset —0.25). 


pendence assuming constant rates of particle production 
and current generation. In contrast, for a small dura¬ 
tion, it is very difficult to predict what should happen a 
priori , and this is why we must perform the numerical 
simulations. 

Our numerical results in Fig. [6] (for N = 8,10,12 cor¬ 
responding to 17 3 ,21 3 ,25 3 ) where B\\ = B± = Eq/2 is 
fixed clearly confirm linear dependence for large dura¬ 
tion (apart from small oscillation which is characteris¬ 
tic to magnetic phenomena). Although j° shows some 
lattice-size dependence (due to the constant shift in the 
ultraviolet region of the produced particle distribution), 
two results for j y with different Ni s are nearly overlaid 
to each other. This clearly indicates that, because the 
constant shift drops off in f d^q, j y has no contamina¬ 
tion from the ultraviolet fluctuations. 

We find for small duration a region of delay in which 
j° and j v remain almost vanishing. The presence of such 
delay is intuitively understandable: it should take some 
time for the wave-function to transform from the free 
plane wave to the Landau-quantized one after a sudden 
switch-on of the magnetic field. We would emphasize 
that what seems non-trivial in Fig. [6] is that the delay 
for j v is more than three times larger than that for j°. 
In our interpretation this difference in the onsets comes 
from the difference in the rates. In Fig. [6] we showed the 
analytically estimated rates (by the dotted line) using 
the formula in Ref. m with a common offset —0.25. 
The onset for j y is more delayed than that for j° since 
its rate dtj v is smaller. 

Finally we briefly comment on the sensitivity when we 
change the box size, i.e., the infrared extrapolation. We 
kept the lattice spacing (and so E 0 a 2 ) and changed the 
lattice size to find that there is a subtle feature in the 


scaling behavior. Since we need to go far into techni¬ 
cal details to address this problem, we will leave it for 
another publication that will be focused on the lattice 
formulation and the systematic checks of the lattice-size 
dependence. 

VI. SUMMARY AND OUTLOOK 

We formulated the production of particles and anti¬ 
particles and checked its validity by looking for a right- 
handed particle excess on CP-odd backgrounds. We lim¬ 
ited our main concern to the CME current at present, 
but this type of the calculation should be applicable for 
a wider range of physical problems. For example, the 
(pair) particle creation at the horizon of the acoustic 
blackhole (see Ref. [323 for a pedagogical article) attracts 
attention, and it would be feasible to combine it with 
CP-odd (or T irreversible) background effects, for which 
we can apply our method. Of course, as mentioned in 
the beginning, possible applications should cover prob¬ 
lems of sphaleron-like transitional processes in the elec- 
troweak and the strong interactions on top of arbitrary 
gauge fields including non-topological ones. 

Here we emphasize the impact of the agreement be¬ 
tween our numerical results and the analytical expres¬ 
sions. It is very non-trivial that a direct evaluation of the 
current according to Eq. § could be consistent with the 
quantum anomaly. Also we emphasize the importance 
of real-time character of the chiral magnetic current gen¬ 
eration. The physical interpretation and the theoretical 
estimate of p -5 are quite problematic and sometimes mis¬ 
leading in the literature. Our successful simulation is the 
first step toward a realistic simulation in experimental 
setups without /iq in the system. 

In this kind of approach to solve the equation of mo¬ 
tion, the anomaly arises from the infrared singularity 
around \p A \ = 0. The major part of the discretization 
error might be attributed to underestimating the singu¬ 
lar contributions with coarse mesh. This implies that 
there may be a hybrid way to extract the singular terms 
analytically and to calculate non-singular parts numeri¬ 
cally m- This would reduce the lattice-size dependence 
and we can make theoretical estimates with improved 
reliability. Though we have confirmed that our CME 
current is close to the one in the continuum limit, it 
would be worth developing such a hybrid algorithm. An¬ 
other direction for the improvement is to include back- 
reaction from gauge fluctuations. Once the back-reaction 
is taken into consideration, it would capture the effect of 
the chiral plasma instability [33] ■ These are works under 
progress 131] . 
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